BCS-BEC crossover at finite temperature in the broken-symmetry phase 
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The BCS-BEC crossover is studied in a systematic way in the broken-symmetry phase between 
zero temperature and the critical temperature. This study bridges two regimes where quantum and 
thermal fluctuations are, respectively, important. The theory is implemented on physical grounds, 
by adopting a fermionic self-energy in the broken-symmetry phase that represents fermions coupled 
to superconducting fluctuations in weak coupling and to bosons described by the Bogoliubov theory 
in strong coupling. This extension of the theory beyond mean field proves important at finite 
temperature, to connect with the results in the normal phase. The order parameter, the chemical 
potential, and the single-particle spectral function are calculated numerically for a wide range of 
coupling and temperature. This enables us to assess the quantitative importance of superconducting 
fiuctuations in the broken-symmetry phase over the whole BCS-BEC crossover. Our results are 
relevant to the possible realizations of this crossover with high-temperature cuprate superconductors 
and with ultracold fermionic atoms in a trap. 

PACS numbers: PACS numbers: 03.75.Ss, 03.75.Hh, 05.30.Jp 



I. INTRODUCTION 



In the BCS to Bose-Einstein condensation (BEG) 
crossoveri242i^i^iSiL2i2ii£, largely overlapping Cooper 
pairs smoothly evolve into non-overlapping composite 
bosons as the fermionic attraction is progressively in- 
creased. These two physical situations (Cooper pairs vs 
composite bosons) correspond to the weak- and strong- 
coupling limits of the theory, while in the interesting 
intermediate-coupling regime neither the fermionic nor 
the bosonic properties are fully realized. Under these cir- 
cumstances, the theory is fully controlled on the weak- 
and strong-coupling sides, while at intermediate coupling 
an interpolation scheme results (as for all crossover ap- 
proaches). These physical ideas are implemented, in 
practice, by allowing for a strong decrease of the chem- 
ical potential at a given temperature when passing from 
the weak- to the strong-coupling limit. 

The BCS-BEC crossover can be considered both be- 
low (broken-symmetry phase) and above (normal phase) 
the superconducting critical temperature. In particu- 
lar, in the normal phase preformed pairs exist in the 
strong-coupling limit up to a temperature T* correspond- 
ing to the breaking of the pairs, while coherence among 
the pairs is established when the temperature is low- 
ered below the superconducting critical temperature Tc- 
This framework could be relevant to the evolution of the 
properties of high-temperature cuprate superconductors 
from the overdoped (weak-coupling) to the underdoped 
(strong-coupling) regions of their phase diagram^ ^. The 
BCS-BEC crossover can be also explicitly realized with 
ultracold fermionic atoms in a trap, by varying their mu- 
tual effective attractive interaction via a Fano-Feshbach 
resonancei^. 

The BCS-BEC crossover has been studied extensively 
in the past, either at T = or for T > T^. At T = 0, the 
solution of the two coupled BCS (mean-field) equations 
for the order parameter A and the chemical potential 



fi has been shown to cross over smoothly from a BCS 
weak-coupling superconductor with largely overlapping 
Cooper pairs to a strong-coupling superconductor where 
tightly-bound pairs are condensed in a Bose-Einstein (co- 
herent) ground statei*2ii2,. For this reason, the BCS mean 
field has often been considered to be a reliable approx- 
imation for studying the whole BCS-BEC crossover at 
T ~ 0. At finite temperature, the increasing importance 
in strong coupling of the thermal excitation of collective 
modes (corresponding to noncondensed bosons) was first 
pointed out by Nozieres and Schmitt-Rink'^ . By their 
approach, the expected result that the superconducting 
critical temperature should approach the Bose-Einstein 
temperature Tbe in strong coupling was obtained (com- 
ing from above Tc) via a (first-order) inclusion of the t- 
matrix self-energy in the fermionic single-particle Green's 
function. The same type of t-matrix approximation (also 
with the inclusion, by some authors, of self-consistency) 
has then been widely adopted to study the BCS-BEC 
crossover above Tc, both for continuum and lattice mod- 
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Despite its conceptual importance, a systematic study 
of the BCS-BEC crossover in the temperature range 
< T < Tc is still lacking. A diagrammatic the- 
ory for the BCS-BEC crossover that extends below Tc 
the self-consistent i-matrix approximation was proposed 
some time ago by Haussmann«. The ensuing coupled 
equations for the order parameter and chemical po- 
tential were, however, solved explicitly only at Tc^ 
leaving therefore unsolved the problem of the study of 
the whole temperature region below Tc- The work by 
Levin and coworkers^S', on the other hand, even though 
based on a "preformed-pair scenario" , has focused mainly 
on the weak-to-intermediate coupling region, where the 
fermionic chemical potential remains inside the single- 
particle band. An extension of the self-consistent t- 
matrix approximation to the superconducting phase for a 
two-dimensional lattice model was considered in Ref. 123. 



In that paper, however, the shift of the chemical poten- 
tial associated with the increasing coupling strength was 
ignored, by keepingft fixed at the noninteracting valueiSi 
The results of Ref. 20 are thus not appropriate to address 
the BCS-BEC crossover, for which the renormalization of 
the chemical potential (that evolves from the Fermi en- 
ergy in weak coupling to half the binding energy of a 
pair in strong coupling) plays a crucial role^'-''^. Addi- 
tional studies have made use of a fermion-boson modelS^, 
especially in the context of trapped Fermi gases^i. 

Purpose of the present paper is to study the BCS-BEC 
crossover in the superconducting phase over the whole 
temperature range from T = to T = Tc, thus filling a 
noticeable gap in the literature. We will consider a three- 
dimensional continuum model, for which the fermionic 
attraction can be modeled by a point-contact interaction. 
As noted in Refs. and [lO, with this model the struc- 
ture of the diagrammatic theory for the single-particle 
fermionic self-energy simplifies considerably, since only 
limited sets of diagrammatic structures survive the reg- 
ularization of the contact potential in terms of the 
fermionic two-body scattering length apmSi^ The di- 
mensionless interaction parameter {kFaF)~^ (where the 
Fermi wave vector kp is related to the density via n = 
kp/^Sn'^)) then ranges from — cxd in weak coupling to +cx) 
in strong coupling. The crossover region of interest is, 
however, restricted in practice by {kF\aF\)~^ ^ 1- 

For this model, a systematic theoretical study of the 
evolution of the single-particle spectral function in the 
normal phase from the BCS to BEC limits has been 
presented recentlyS^. Like in Ref. y, also in Ref. |2^ 
the coupling of a fermionic single-particle excitation to 
a (bosonic) superconducting fluctuation mode was taken 
into account by the i-matrix self-energy. This approx- 
imation embodies the physics of a dilute Fermi gas in 
the weak-coupling limit and reduces to a description 
of independent composite bosons in the strong-coupling 
limit. In this way, single-particle spectra were obtained 
in Ref. |2^ as functions of coupling strength and temper- 
ature. 

In the present paper, the t-matrix approximation for 
the self-energy is suitably extended below Tc. In par- 
ticular, the same superconducting fluctuations, that in 
Refs. y and |2^ were coupled to fermionic independent- 
particle excitations above T^, are now coupled to 
fermionic BCS-like single-particle excitations below Tc- 
In the strong-coupling limit, it turns out that these super- 
conducting fluctuations merge in a nontrivial way^^ into 
a state of condensed composite bosons described by the 
Bogoliubov theory, and evolve consistently into a state of 
independent composite bosons above Tc (as the Bogoli- 
ubov theory for point-like bosons doesSi). In this way, a 
direct connection is established between the structures of 
the single-particle fermionic self-energy above and below 
Tci as they embody the same kind of bosonic mode which 
itself evolves with temperature. 

A comment on the validity of the Bogoliubov theory at 
finite temperature (and, in particular, close to the Bose- 



Einstein transition temperature Tbe) might be relevant 
at this point. A consistent theory for a dilute condensed 
Bose gas was developed long ago in terms of a (small) 
gas parameteriSSiSS, of which the Bogoliubov theory'^° is 
only an approximate form valid at low enough temper- 
atures (compared with Tbe)- That theory correctly de- 
scribes also the dilute Bose gas in the normal phasa^S, 
whereas the Bogoliubov theory (when extrapolated above 
the critical temperature) recovers the independent-boson 
form (albeit in a non-monotonic way, with a discontinu- 
ous jump affecting the bosonic condensate^^). It would 
therefore be desirable to identify (at least in principle) 
a fermionic theory that, in the strong-coupling limit of 
the fermionic attraction, maps onto a more sophisti- 
cated bosonic theory, overcoming the apparent limita- 
tions of the Bogoliubov theory. In practice, however, it 
should be considered already a nontrivial achievement 
of the present approach the fact that the bosonic Bogoli- 
ubov approximation can be reproduced from an originally 
fermionic theory. For these reasons, and also because it 
is actually the intermediate-coupling (crossover) region 
to be of most physical interest, in the following we shall 
consider the Bogoliubov approximation as a reasonable 
limiting form of our fermionic theory. 

As it is always the case for the BCS-BEC crossover 
approach, implementation of the theory developed in the 
present paper rests on solving two coupled equations 
for the order parameter A and the chemical potential 
/i. The equations here considered for A and /z general- 
ize the usual equations already considered at the mean- 
field leve l-^i^i^ , by including fluctuation corrections. Our 
equations reproduce the expected physics in the strong- 
coupling limit, at least at the level of approximation here 
considered. Their solution provides us with the values 
of A and /i as functions of coupling strength {kpaF)~^ 
and temperature T, thus extending results obtained pre- 
viously at the mean-field level. In particular, the order 
parameter is now found to vanish at a temperature (close 
to) Tc even in the strong-coupling limit, while it would 
had vanished close to T* at the mean-field level^. 

The analytic continuation of the fermionic self-energy 
to the real frequency axis is further performed to ob- 
tain the single-particle spectral function A(k, w), that we 
study in a systematic way as a function of wave vector k, 
frequency w, coupling strength {kpCiF)'^, and temper- 
ature T. In this context, two novel sum rules (specific 
to the broken-symmetry phase) are obtained, which pro- 
vide compelling checks on the numerical calculations. In 
addition, the numerical calculations are tested against 
analytic (or semi-analytic) approximations obtained in 
the strong-coupling limit. The study of a dynamical 
quantity like A{\<i,uj) enables us to attempt a compari- 
son with the experimental ARPES and tunneling spectra 
for cuprate superconductors below Tc, for which a large 
amount of data exists showing peculiar features for differ- 
ent doping levels and temperatures. As in Ref. 123 above 
Tc, this comparison concerns especially the experimental 
data about the M points in the Brillouin zone of cuprates, 



where pairing effects are supposed to be stronger than 
along the nodal lines. 

Our main results are the following. About thermo- 
dynamic quantities, we will show that fluctuation cor- 
rections over and above mean field are especially impor- 
tant at finite temperature T "^Tc when approaching the 
strong-coupling limit. At zero temperature, fluctuation 
corrections to thermodynamic quantities turn out to be of 
some relevance only in the intermediate-coupling region. 
This supports the expectation^ that the BCS mean field 
at zero temperature should describe rather well the BCS- 
BEC crossover essentially for all couplings. Regarding in- 
stead dynamical quantities like A{\<i,oj), our calculation 
based on a "preformed-pair scenario" reveals two distinct 
spectral features for w < 0. These features, which have 
different temperature and doping dependences, together 
give rise to a peak-dip-hump structure which is actively 
debated for the ARPES spectra of cuprate superconduc- 
tors. Our results differ from those previously obtained by 
other calculationsi^ also based on a "preformed-pair sce- 
nario" , where a single feature was instead obtained in the 
spectral function for w < 0. An explanation of this dis- 
crepancy between the two calculations will be provided. 
It will also turn out from our calculation that the co- 
herent part of v4(k, to) for a; < follows essentially a 
BCS-like behavior as far as its wave-vector dependence 
is concerned, albeit with a gap value which contains an 
important contribution from fluctuations at finite tem- 
perature. The same BCS-like behavior is not found, 
however, by our calculation for the dependence of the 
spectral weight of the coherent peak on temperature and 
coupling. This evidences a dichotomy in the behavior of 
v4(k, w), according to which of its dependences one is af- 
ter. Such a dichotomy is clearly observed in experiments 
on cuprate superconductors, in good qualitative agree- 
ment with the results obtained by our calculations.'^^. A 
detailed quantitative comparison of our results with the 
experimental data on cuprates would, however, require a 
more refined theoretical model, as to include the quasi- 
two-dimensional lattice structure, the d-wave character of 
the superconducting gap, and also a fermionic attraction 
that depends effectively on doping (and possibly on tem- 
perature). Future work on this subject should address 
these additional issues. 

The present theory could be improved in several ways. 
In the present approach, the effective interaction between 
the composite bosons is treated within the Born approx- 
imation. For a dilute system of composite bosons one 
knows how to improve on this result, as shown in Ref.llOl 
(see also Ref. 33). In addition, the Bogoliubov descrip- 
tion for the composite bosons could be also improved, 
for instance, by extending to the composite bosons the 
Popov treatment for point-like bosons^^. Finally, on 
the weak-coupling side of the crossover the BCS theory 
could be modified by including the contributions shown 
by Gor'kov and Melik-BarkhudaroAi^i to yield a finite 
renormalization of the critical temperature and of the 
gap function even in the extreme weak-coupling limit. 



Work along these lines is in progress. 

The plan of the paper is as follows. In Sec. II we discuss 
our choice for the fermionic self-energy in the supercon- 
ducting phase, from which the order parameter A and the 
chemical potential /i are obtained as functions of tem- 
perature and coupling strength, and the spectral func- 
tion A(k, uj) also results. Analytic results are presented 
in the strong-coupling limit, where the order parameter 
is shown to be connected with the bosonic condensate 
density of the Bogoliubov theory. In addition, the an- 
alytic continuation of our expressions for the fermionic 
self-energy and spectral function is carried out in detail. 
In Sec. Ill we present our numerical calculations, and 
discuss the results for the single-particle spectral func- 
tion in the context of the available experimental data for 
high-temperature cuprate superconductors. Section IV 
gives our conclusions. In Appendix A two sum rules are 
derived for the superconducting phase, which are used as 
checks of the numerical results. 



II. DIAGRAMMATIC THEORY FOR THE 

BCS-BEC CROSSOVER IN THE 

SUPERCONDUCTING PHASE 

In this section, we discuss the choice of the fermionic 
single-particle self-energy in the superconducting phase 
for a (three-dimensional) continuum system of fermions 
mutually interacting via an attractive point-contact po- 
tential, with an s-wave order parameter. We shall place 
special emphasis to the strong-coupling limit of the the- 
ory, where composite bosons forms as bound fermion 
pairs. We extend in this way below T^ an analogous treat- 
ment for the self-energy, made previously in the normal 
phase to calculate the single-particle spectral function. ■^^ 

Knowledge of the detailed form of the attractive inter- 
action is not generally required when studying the BCS- 
BEC crossover. Accordingly, one may consider the sim- 
ple form vq5{y) of a "contact" potential, where vq is a 
negative constant. This choice entails a suitable regular- 
ization in terms, e.g., of a cutoff /cq in wave- vector space. 
In three dimensions, this is achieved via the scattering 
length ap of the associated fermionic two-body problem, 
by choosing vq as foUowsiSi: 



Vo = 



mfcn 



map kg 

m being the fermion mass. With this choice, the classifi- 
cation of the (fermionic) many-body diagrams is consid- 
erably simplified not only in the normal phaseiS but also 
in the broken-symmetry phase^^, since only specific dia- 
grammatic substructures survive when the limit fcg — > oo 
(and thus f o ^ 0) is eventually taken. 

In particular, the particle-particle ladder depicted in 
Fig. 1 (a) survives the regularization of the potential J^ It 
is obtained by the matrix inversion: 



rii(g) ri2(g) 



xii(-'7) xi2{q) 
xi2{q) xii(9) 



X [xii(g)xii(-'7)-xi2(g)']-\2) 



with the notation 



Xii(9) 



Awap 



dp 



^Xl^ii(^ + '?)^ii(~^) 



m 



Xi2{q) 



y^2^^E^i2(p + g)e2i(-p) 



(3) 

(4) 



In these expressions, q = (q, ^u) and p = (p, ujn), where 
q and p are wave vectors, and il,y = l-nvj (i {v in- 
teger) and w„ = (2n + l)7r//3 {n integer) are bosonic 
and fermionic Matsubara frequencies, respectively (with 
/3 = {kBT)~^, ks being the Bohzmann's constant); 



gil{p,UJn) = - 
^2l(P,W„) 



e(p) 



E{pr 

A 



Eip) 



-ujt, 



= -^22(-p, 
^12(P,W„) 



(5) 



are the BCS single-particle Green's functions in Nambu 
notation, with ^(p) = p^/(2m) — /i and E{p) = 
\/C(p)^ + ^^ for an isotropic (s-wave) order parameter 
A. [Hereafter, we shall take the order parameter to be 
real with no loss of generality.] 

The expressions ^ and Q for Xii((z) and Xi2(9) con- 
siderably simplify in the strong- coupling limit (that is, 
when /3/i — > — oo and A ^ |/i|). In this limit, one then 
obtains for the matrix elements (j21>^*^: 



TiiCg) = r22(-g) 



Stt /^b -t- zf}^ + qV(4m) 



rn^ap 



Esicif - m,f 



and 

ri2(g) 

where 



r2i(<z) 



Stt 



A^B 



m2aF £^i3(q)2 - {i^uf 



Eeici) = 




\2mB 



M-B 



mI 



(6) 

(7) 

(8) 



has the form of the Bogoliubov dispersion relation^° 
{mB — 2m being the bosonic mass, ps = ^^/(4|a*|) = 
2/^-f eo the bosonic chemical potential, and eg = {ma^p)~^ 
the bound-state energy of the associated fermionic 
two-body problem). The above relation between the 
fermionic and bosonic chemical potentials holds provided 
A^B ^ eo (cf. also Sec. IID). Note that ^mb can be cast in 
the Bogoliubov form 



Mb = ^2(0) rio(r) 



(9) 



where ^2(0) = A^nap/Tn is the residual bosonic 
interaction^Jfi and no(r) = /:^^{T)m'^ap/{Si:) is the con- 
densate density. The relation (Q is formally obtained al- 
ready at the (BCS) mean-field level^^, albeit with an un- 
specified dependence of no(T') on temperature. Within 



q-k q-k' q-k q-k' q-k q-k' 



'L^ 



'l^ 



r(q) 



■»-'r 'l- 



(a) 



k' 



k k' 



'r h > » ) « ) (r 



k k' 

'l 'r 



(b) 



'l^ 



r(q) 



(c) 




k' (d) 



FIG. 1: (a) Particle-particle ladder in the broken-symmetry 
phase. Conventions for four-momenta and Nambu indices are 
specified. Dots delimiting the potential (broken line) repre- 
sent ra Pauli matrices. Only combinations with II = £'l and 
£r — I'll occur owing to the regularization we have adopted for 
the potential, (b) Fermionic self-energy diagram associated 
with the expression 1121 in the normal phase, (c) Fermionic 
self-energy diagram associated with the expressions 11311 and 
11411 in the broken-symmetry phase, (d) BCS contribution 
(15) to the self-energy. 



our fiuctuation theory, the temperature dependence of 
no (r) will coincide in strong coupling with the expression 
given by the Bogoliubov theory (see Sec. IID). In partic- 
ular, at zero temperature and at the lowest order in the 
residual bosonic interaction^^, no reduces to the bosonic 
density ub = n/2 and fiB is given by 2k'pap j {iirm) . 

Note further that the above result for V2 (0) can be cast 
in the bosonic form ^2(0) — AnaB/iTiB with ub ~ 2a p. 
The present theory thus describes the effective interac- 
tion between the composite bosons within the Born ap- 
proximation, while improved theoriesiS*^ for qb would 
give smaller values for the ratio aB/ap- These improve- 
ments will not be considered in the present paper. 

Apart from the overall factor —8Tr/{m?ap) (and a sign 
difference in the off-diagonal componen1i2&), the expres- 
sions (0 and 10 coincide with the normal and anoma- 
lous non-condensate bosonic Green's functions within the 
Bogoliubov approximation^^, respectively. These expres- 
sions will be specifically exploited in Sec. IID, where the 
strong-coupling limit of the fermionic self-energy will be 
analyzed in detail. 

In the normal phase, on the other hand, the BCS 
single-particle Green's functions are replaced by the bare 
single-particle propagator ^o(p) = [i'^n — £,ip)]^^i while 
for arbitrary coupling the particle-particle ladder ac- 
quires the form: 



ro('7) = - 



- /. 



\ 4nap 



dp 

(2^ 



tanh(/3^(p)/2) + tanh(/3^(p - q)/2) m 

2(C(p) + e(p-g)-*f^.) P^ 



(10) 



In particular, in the strong-coupling limit the expression 
pO|l reduces to 



ro(<?) 



Stt 



1 



m?aF i^u — q^/(4TO) 



(11) 



which coincides (apart again from the overall factor 
— 87r/(m^ai?)) with the free-boson Green's function. 

The above quantities constitute the essential ingre- 
dients of our theory for the fermionic self-energy and 
related quantities in the broken-symmetry phase. As 
shown in Ref. 26| they also serve to establish a mapping 
between the fermionic and bosonic diagrammatic struc- 
tures in the broken-symmetry phase, in a similar fashion 
to what was done in the normal phase^fi. 



A. Choice of the self-energy 

In a recent study^^ of the single-particle spectral 
function in the normal phase based on the BCS-BEC 
crossover approach, the fermionic self-energy was taken 
of the form: 

So(fc) = -^Ero('?)^o(g-fc) (12) 

q 

where V is the quantization volume and k = (k, lj^) 
is again a four-vector notation with wave vector k and 
fermionic Matsubara frequency uog (s integer). In this ex- 
pression, Fq ((?) is given by Eq. H10() for arbitrary coupling 
and Go{k) is the bare single-particle propagator. The self- 
energy diagram corresponding to the expression H12|) is 
depicted in Fig. 1(b). The fermionic single-particle ex- 
citations are effectively coupled to a (bosonic) supercon- 
ducting fluctuation mode, which reduces to a free com- 
posite boson in the strong-coupling limit. Physically, the 
choice (I12II for the self-energy entails the presence of a 
pairing interaction above Tc, which can have significant 
influence on the single-particle (as well as other) proper- 
ties. 

In the present paper, we choose the self-energy in the 
broken-symmetry phase below Tc^ with the aim of re- 
covering the expression H12|l when approaching Tc from 
below and the Bogoliubov approximation for the compos- 
ite bosons in the strong-coupling limit. To this end, we 
adopt the simplest approximations to describe fermionic 
as well as bosonic excitations in the broken-symmetry 
phase, which reduce to bare fermionic and free bosonic 
excitations in the normal phase, respectively. These 
are the BCS single-particle Green's functions © (in the 
place of the bare single-particle propagator Qq) and the 
particle-particle ladder (|2J) (in the place of its normal- 
phase counterpart Fq). By this token, the fermionic self- 
energy H12|l is replaced by the following 2x2 matrix: 

I]fi(A;)--E^2(-fc) = -^^rn(g)t/n(g-fc) (13) 



^Uk)^^Uk)^-^Y.^i2{q)gi2{q'k) (14) 

where the label L refers to the particle-particle ladder. 
The corresponding self-energy diagram is depicted in 
Fig. l(c)^ 

The choice (|13|) and (|14|l for the self energy is made 
on physical grounds. A formal "ab initio" derivation of 
these expressions can also be done in terms of "conserv- 
ing approximations" in the Baym-Kadanoff sense, that 
hold even in the broken-symmetry phase^^. In such a 
formal derivation, however, the single-particle Green's 
functions entering Eas. l|13|l and (|14|l (also through the 
particle-particle ladder jSJ) would be required to be self- 
consistently determined with the sam,e self-energy inser- 
tions. In our approach, we take instead the single-particle 
Green's functions to be of the BCS form 10 • The or- 
der parameter A and chemical potential /i are obtained, 
however, via two coupled equations (to be discussed in 
Sec. lie) that include the self-energy insertions \i'A\ and 
(|14|l . In this way, we will recover the Bogoliubov form 
© and (UI) for the particle-particle ladder not only at 
zero temperature but also at flnite temperatures (and, in 
particular, close to the Bose-Einstein transition temper- 
ature) . 

The choice (|13|l and (|14() for the self energy is not ex- 
haustive. In the broken-symmetry phase there, in fact, 
exists an additional self-energy contribution that survives 
the regularization (^) of the interaction potential in the 
limit /cq ~* OO; even though it does not contain particle- 
particle rungs'^^. This additional self-energy diagram is 
the ordinary BCS contribution depicted in Fig. 1(d), with 
the associated expression 



BCS( 



Sf2^-(A:) = I]2TnA^) = -A 



(15) 



while the corresponding (Hartree-Fock) diagonal ele- 
ments vanish with the regularization we have adopted. 
Relating the expression ((T3|) to the diagram of Fig. 1(d) 
rests on the validity of the BCS gap equation [Eq. H2U|) 
below] , for arbitrary values of the chemical potential. For 
this, as well as for an additional reason (cf. Sec. IID), we 
shall consistently consider that equation to hold for the 
order parameter A. 

The choice ((TK)l alone would be appropriate to describe 
the system in the weak-coupling (BCS) limit, where the 
superconducting fluctuation contributions \V6\ and l|14l) 
represent only small corrections. In the intermediate- 
and strong-coupling regions, on the other hand, both 
contributions (|13|I - H14I) and IjlSII might become equally 
significant (depending on the temperature range below 
Tc). We thus consider both contributions simultaneously 
and write the fermionic self-energy in the matrix form: 



Sii(A:) 
^2i{k) 



Sl2(fc) 
S22(fc) 

Sfi(A:) 



Sfi(fc) 



'^BCS 



{k) 



^Uk) 



^BCS 
■^12 



^221'^) 



(fc) 



(16) 



In the following, however, we shall neglect Sf'2 in com- 
parison to Y.f^^ . It will, in fact, be proved in Sec. IID 
that, in strong coupling, Efj is subleading with respect 
to both 'Sf'f^ and 'S^i- Inclusion of S^j is thus not re- 
quired to properly recover the Bogoliubov description for 
the composite bosons in the strong-coupling limit. 

To summarize, the fermionic single-particle Green's 
functions are obtained in terms of the bare single-particle 
propagator Go{k) and of the self-energy l(T^ and l(TB|) via 
the Dyson's equation in matrix form: 



G21 (k) G22 (k) 



Goik)- 






-Go{-ky 



(fc) 



vscs 



'^22ik) 



(17) 



If only the BCS contribution H15() to the self-energy were 
retained, the fermionic single-particle Green's functions 
Gij{k) {i,j = 1,2) would reduce to the BCS form JSJ. 
Upon including, in addition, the fluctuation contribution 
p3|l to the self-energy, modified single-particle Green's 
functions result, which we are going to study as functions 
of coupling strength and temperature. 



B. Comparison with the Popov approximation for 
dilute superfluid fermions 

The choice of the self-energy (|13|l and H15() resem- 
bles the approximation for the self-energy introduced 
by Popov'^^ for superfluid fermions in the dilute limit 
A;F|a_F| ^ 1 (with ap < 0). There is, however, an im- 
portant difference between the Popov fermionic approx- 
imation and our theory. We include in Eq. (|13|l the full 
Fii obtained by the matrix inversion of Eq. (|2Jl; Popov in- 
stead neglects X12 therein and approximate Fn by l/xn, 
thus removing the feedback of the Bogoliubov- Anderson 
mode on the diagonal fermionic self-energy Sn. Retain- 
ing this mode is essential when dealing with the BCS- 
BEC crossover, to describe the composite bosons in the 
strong-coupling limit by the Bogoliubov approximation, 
as discussed in Sec. IIA. Approaching the weak-coupling 
limit, on the other hand, the presence of the Bogoliubov- 
Anderson mode becomes progressively irrelevant and the 
self-energies coincide in the two theories. As a check on 
this point, we have verified that, in the weak-coupling 
limit and at zero temperature, En obtained by our the- 
ory (using the numerical procedures discussed in Sec. Ill) 
reduces to 4:napn/{2m), which is the expression obtained 
also with the Popov approximation^>2, in the absence of 
the Bogoliubov- Anderson mode. 

There is another difference between the Popov 
fermionic approximation and our theory as formulated in 
Sec. IIA, which concerns the off-diagonal fermionic self- 
energy S12. Our expression (|15|l for E12 was obtained 
from the diagram of Fig. 1(d), where the single particle 
line represents the off-diagonal BCS Green's function of 
Eq. (O with no insertion of the diagonal self-energy En. 



Within the Popov approximation, on the other hand, E12 
is defined formally by the same diagram of Fig. 1(d), 
but with the single-particle line being fully self-consistent 
(and thus including En). Since En turns out to ap- 
proach a constant value Eq in the weak-coupling limit 
(as discussed above), inclusion of En ~ Eq can be sim- 
ply made by a shift of the chemical potential (such that 
fi —f ^ — Eg). This shift affects, however, the value of the 
gap function A in a non-negligible way even in the ex- 
treme weak-coupling limit. Neglecting this shift, in fact, 
results in a reduction by a factor e^^^ of the BCS asymp- 
totic expression (8ei;-/e^) exp[7r/(2fci?aF)] for A (where 
ep — kp/{2m)). Inclusion of the shift Eo is thus im- 
portant to recover the BCS value for A in the (extreme) 
weak-coupling limit. 

The need to include the constant shift Eg on the weak- 
coupling side of the crossover was also discussed in Ref.l25l 
while studying the spectral function A(k, uj) in the nor- 
mal phase with the inclusion of pairing fluctuations. In 
that context, inclusion of the shift Eq proved necessary to 
have the pseudogap depression of A(k, uj) centered about 
UJ = 0. Inclusion of the shift Eq in the broken-symmetry 
phase (at least when approaching the critical tempera- 
ture from below) is thus also necessary to connect the 
spectral function A(k, oj) with continuity in the weak- 
coupling side of the crossover. 

Combining the above needs for A and A(k, cu), we have 
introduced the constant shift Eq for all temperatures be- 
low Tc, by replacing /i with /i — Eq in the BCS Green's 
functions (0) entering the convolutions Q and Q). The 
same replacement is made in the gap equation [Eq. H2()|l 
below]. In the Dyson's equation H17|l . however, fj, is left 
unchanged since the constant shift Eg is already con- 
tained in Ell (A;) as soon as its fc-dependence is irrelevant. 
Accordingly, we have included this constant shift in the 
calculation of both thermodynamic and dynamical quan- 



tities in the weak-coupling side for [kpap)' 



< 



-0.5, 



and neglected it for larger couplings when Eii(fc) can no 
longer be approximated by a constant. 

It turns out that the temperature dependence of Eg is 
rather weak in the above coupling range. A plot of Eg vs 
T/Tc and {kpap)~^ is shown in Fig. [3 Here, the criti- 
cal temperature Tc is obtained by applying the Thouless 
criterion from the normal phase as was done in Ref. |25| 
(this procedure to obtain Tc will be used in the rest of 
the paper). In this plot, the constant shift Eg is obtained 
as Eg = ReEn(|k| = yj2m,{\i — Eg), w = 0), in analogy 
to what was also done in Ref. l25[ Here, En(k, cj) is the 
analytic continuation to the real frequency axis of the 
Matsubara self-energy Eii(k, Wg) discussed in Sec. HE. 



C. Coupled equations for the order parameter and 
the chemical potential 

Thermodynamic quantities, such as the order param- 
eter A and the chemical potential ^, are obtained di- 
rectly in terms of the Matsubara single-particle Green's 




(kpap) 



equation has the same formal structure of the BCS gap 
equation, although the numerical values of the chemi- 
cal potential entering Eq. (|2()|l differ from those obtained 
by the BCS density equation. This procedure ensures 
that the bosonic propagators (0) in the broken-symmetry 
phase are gapless, as shown explicitly by the Bogoliubov- 
type expressions (0 and Q in the strong-coupling limit. 
In general, in fact, there is no a priori guarantee that 
a given (conserving) approximation for fermions would 
result into a "gapless" approximation^^ for the compos- 
ite bosons in the strong-coupling limit of the fermionic 
attraction. 

Including fluctuation corrections to the BCS density 
equation as in Eq. ifT^ . on the other hand, results in 
the emergence of important effects in the strong-coupling 
limit of the theory, as discussed next. 



FIG. 2 

T (units of T, 



Self-energy shift Eq (in units of ep) vs temperature 
of Tc) and coupling {kpap)'^ ■ 



functions, without the need of resorting to the analytic 
continuation to the real frequency axis. 

Quite generally, the order parameter A is defined in 
terms of the "anomalous" Green's function Gi2(k, Wj) via 
A = vo{'ip^{r)ipi{r)) [cf. Eq. HA12p ]. where the strength 
Wo of the contact potential is kept to comply with a stan- 
dard definition of BCS theory^*^. One obtains: 



A 



■Vo 



jj^-pT.Gi2{K 



^s) 



(18) 



By the same token, the chemical potential /i can be 
obtained in terms of the "normal" Green's function 
Gii(k, Ws) via the particle density n: 



n = 2 



dk 1 



(27r)3 (3 






'Gii(k,c^,) (19) 



where 77 = 0+. The two equations H18|l and H19|l are 
coupled, since the Green's functions depend on both A 
and fj,. The results of their numerical solution will be 
presented in the next section for various temperatures 
and couplings. 

In the following treatment, we shall deal with the two 
equations (|18|l and H19|l on a different footing. Specifi- 
cally, we will enter in the density equation (|19|l the ex- 
pression for the normal Green's function obtained from 
Eq. H17|) . that includes both BCS and fluctuation contri- 
butions (see Eq. H35(l below). We will use instead in the 
gap equation H18|l the BCS anomalous function (jSJ, that 
includes only the BCS self-energy H15|l . In this way, the 
gap equation ((TH|l reduces to the form 



Anai 



dk 



tanh(/3£;(k)/2) m 
2E(k) k2 



(20) 



where the regularization of the contact potential in terms 
of the scattering length ap has been introduced. This 



D. Analytic results in the strong-coupling limit 

We proceed to show that the original fermionic theory, 
as defined by the Dyson's equation (|17() . maps onto the 
Bogoliubov theory for the composite bosons which form 
as bound-fermion pairs in the strong-coupling limit. To 
this end, we shall exploit the conditions /3|/i| 3> 1 and 
A <C \fJ-\ {p < 0) (which define the strong-coupling limit) 
in the (Matsubara) expressions (|13|l and H14|) for E^(p), 
thus also verifying that Ef^ can be neglected. 

These expressions are calculated by performing the 
wave- vector and frequency convolutions with the approx- 
imate expressions (0) and lO for the particle-particle lad- 
der and the expressions jHl for the BCS single-particle 
Green's functions. 

Upon neglecting contributions that are subleading un- 
der the above conditions, we obtain in this way for the 
diagonal part of the self-energy: 



^lliKiOs) 



rn^ap 



dq 

f2^ 



uliq)b{EB{q)) 



lUJs 



vliq)bi-EB{q)) 



■E{q-k)+EBiq) 



Eiq~k)~EB{q) 



(21) 



In this expression, E{k) is the BCS dispersion of Eqs. IjSJl, 
EB{q) is the Bogoliubov dispersion relation (jSJ, b{x) — 
[exp{(3x) — 1]~^ is the Bose distribution, and 



«l(q) = "l(q)-i 



q 

2'mE 



fJ-B - EBiq) 



2EB{q) 



(22) 



are the standard bosonic factors of the Bogoliubov 
transformation'^'' . 

In the numerators of the expressions within brackets in 
Eq. (|21|l . the Bose functions are peaked at about q = 
and vary over a scale q^/(2ms) ~ T ^ |/i|. Similarly, 
the factors M^(q) and Vg{q) are also peaked at about 
q = and vary over a scale q^/{2mB) ~ A's <C \fJ,\. The 
denominators in the expression H21|) . on the other hand. 



vary over the much larger scale \fj,\. For these reasons, we 
can further approximate the expression 1)21(1 as follows: 



where 



n'siT) 



Sfi(k,^,) 



dq 



Stt 
mPap 



tOJs 



■?(k) 



,(T) (23) 



(27r)= 



[uliq)biEB{<D)~vUcL)bi-EBi<D)] 

(24) 

identifies the bosonic noncondensate density according 
to Bogoliubov theorji^. Note that in the normal phase 
(when the condensate density no{T) of Eq. Q vanishes), 
the noncondensate density (|24|l becomes the full bosonic 
density ub — n/2 and Eq. (|^ reduces to the expression 
obtained in Ref. [13 directly from the form (|12|l of the 
fermionic self-energy. 

The off-diagonal self-energy S^2(^) can be analyzed in 
a similar way. Since its magnitude is supposed to be the 
largest at zero temperature, we estimate it correspond- 
ingly for k = and w^ = as follows: 



Ef2(k = 0,L^, = 0) 



Stt ^iB^ 
m^aF 2 
1 



dk' 



1 



(2^)3S(k')^i3(k') 



{E{k') + EB{k')) 



(25) 



At the leading order, we can neglect both A and /is in 
the integrand, where the energy scale |/i| dominates. We 
thus obtain Sf'2(fc = 0) « A(A2/(2|^p)) in the strong- 
coupling limit (where the relation /i^ = A^/(4|/i|) - see 
below - has been used). This represents a subleading 
contribution in the small dimensionless parameter A/|/z| 
with respect to both the BCS contribution S^2^^(fc) = 
—A and the diagonal fluctuation contribution Ef]^(fc). It 
can accordingly be neglected. 

Within the above approximations, the inverse p7(l of 
the fermionic single-particle Green's function reduces to: 

G^iik) G^,\k) 



e(k) 



s+a^) 



A 



lUJs 



^(k) 



<W) 



with the notation 



^0 — 



Stt 
m?aF 



n'siT) . 



1(26) 



(27) 



From Eq. H2()|l we get the desired expression for 
G'ii(k, Ws) in the strong-coupling limit: 



Gll(\i,LL!s) 






(28) 



where we have discarded a term of order Aq/|^| with 
respect to |^|. Note that Eq. (|28|) has the same formal 
structure of the corresponding BCS expression |SJl, with 



the replacement £;(k) -^ E(k) ^ ^^(k)^ + (A2 + A^ 



We rewrite it accordingly as: 

„ .. . ^^(k) 

Gii{k, ujs) ~ 



^(k) 



-Eik) 



E{k) 



(29) 



with the modified BCS coherence factors ^^(k) = 1 — 
^2(k) = (l-C(k)/^(k))/2. 

Before making use of the asymptotic expression H29() 
in the density equation p9|l . it is convenient to manipu- 
late suitably the gap equation ((2()() in the strong-coupling 
limit. Expanding 1/-B(k) therein as [1- AV(2f (k)]/^(k) 
and evaluating the resulting elementary integrals, one ob- 
tains: 



^ c. 2 (yn^ - 2 |a*i 



(30) 



Setting further 2/i = — eo + /is, one gets the relation 
A2/(4|^|) = fiB quoted already after Eqs. (jHJ and ^^. 
Let's now consider the density equation H19|l . With the 
BCS-like form H29|) one obtains immediately: 



dk 

^2^ 



f}^(k) 



(31) 



that holds for T ^ eo, at temperatures well below the 
dissociation threshold of the composite bosons. Similarly 
to what was done to get the gap equation ^U^, in Eq. J^ 
one expands l/£;(k) as [1 - (A^ -|- Ag)/(2C(k)]/C(k) and 
evaluates the resulting elementary integrals, to obtain: 



^ ^P /a2 , a2 



An 



(A' + Ag) 



(32) 



Recalling the definition (|27|l for Aq, as well as the ex- 
pressions (|30|l and for the order parameter, which we 
rewrite in the form 



A^ = 4^ no(T) 



(33) 



in analogy to Eq. ((77|) . the result (|S2)l becomes eventu- 
ally: 



n = 2(n's(T) + no(T)) 



(34) 



that holds asymptotically for T ^ eo. 

These results imply that, in the strong-coupling limit, 
the original fermionic theory recovers the Bogoliubov 
theory for the composite bosons, not only at zero tem- 
perature but also at any temperature in the broken- 
symmetry phase. Accordingly, the noncondensate den- 
sity n'g{T) is given by the expression H24() . the bosonic 
factors Vg{q) and u^(q) are given by Eq. H22|l. and the 
dispersion relation -E_B(q) is given by Eq. (jHJ. In the 
strong-coupling limit, the present fermionic theory thus 
inherits all virtues and shortcomings of the Bogoliubov 
theory for a weakly-interacting Bose gasSI. The present 
fermionic theory at arbitrary coupling then provides an 
interpolation procedure between the Bogoliubov theory 



for the composite bosons and the weak-couphng BCS the- 
ory plus pah'ing fluctuations. Both these analytic limits 
will constitute important checks on the numerical calcu- 
lations reported in Sec. III. Note that inclusion of the 
off-diagonal fluctuation contribution Sf'2(^) to the self- 
energy is not required to recover the Bogoliubov theory 
in strong coupling. For this reason, we will not consider 
I]i2{k) altogether in the numerical calculations presented 
in Sec. Ill, as anticipated in Eq. (|17|l . 

The above analytic results enable us to infer the 
main features of the temperature dependence of the or- 
der parameter in the strong-coupling limit. In partic- 
ular, the low-temperature behavior r7,o(T) = 7^o(0) ~ 
m_B(fc_BT)^/(12c)( where c = \/noV2{0)/mB is the sound 
velocity) within the Bogoliubov approximation, implies 
that A(T) decreases from A(0) with a T^ behavior, in 
the place of the exponential behavior obtained within the 
BCS theory (with an s-wave order parameter)^. In addi- 
tion, in the present theory the order parameter vanishes 
over the scale of the Bose-Einstein transition tempera- 
ture Tbe, while in the BCS theory it would vanish over 
the scale of the bound-state energy eo of the composite 
bosons. 

Note finally that the fermionic quasi-particle disper- 
sion -E(k), entering the expression (|29|l of the diagonal 
Green's function in the strong-coupling limit, contains 
the sum A^ -I- Aq instead of the single term A^ of the 
BCS dispersion £;(k). 



given coupling from matrix inversion of Eq. H17|l : 



Gii(k,Ws) 



iuJs -C(k) - i;ii(k,Ws) 
A2 



X (icos + ^(k) - S22(k, Us)) ~ A2]-i 



(35) 



(36) 



Consider first the normal Green's function (|35|l . which 
we rewrite in the compact form 



Gii(k,^s) 



1 



iujs -^(k) -CTii(k,Ws) 
with the short-hand notation 

A2 



(37) 



C7ll{k,UJs) = Ell(k,LJs) 



e(k)- I]22(k,C^. 



(38) 

To perform the analytic continuation of this expression, 
we look for a function crii(k, z) of the complex frequency 
z which satisfies the following requirements at any given 
k: 

(i) It is analytic off the real axis; 

(ii) It reduces to CTii(k, Ws) given by Eq. lp!H|l when z 
takes the discrete values icog on the imaginary axis; 
(iii) Its imaginary part is negative (positive) for Im z > 
(Imz < 0) ; 

(iv) It vanishes when |z| — > oo along any straight line 
parallel to the real axis with Im z ^ 0. 

Once the function i7ii(k, z) is obtained, the expression 



E. Spectral function and sum rules 

We pass now to identify the form of the spectral func- 
tion A(k, u) associated with the approximate choice of 
the Matsubara self-energy of Eq. (|17|1 . To this end, we 
need to perform the analytic continuation in the complex 
frequency plane, thus determining the retarded fermionic 
single-particle Green's functions from their Matsubara 
counterparts. The approach developed in this subsec- 
tion holds specifically for the approximate choice for the 
self-energy of Eq. H17|l . It thus differs from the general 
analysis presented in the Appendix which holds for the 
exact Green's functions, irrespective of any specific ap- 
proximation. 

In general, the process of analytic continuation to 
the real frequency axis from the numerical Matsub- 
ara Green's functions proves altogether nontrivial, as 
it requires in practice recourse to approximate nu- 
merical methods such as, e.g., the method of Fade 
approximants^. We then prefer to rely on a proce- 
dure whereby the analytic continuation to the real fre- 
quency axis is achieved by avoiding numerical extrapola- 
tions from the Matsubara Green's functions. 

The fermionic normal and anomalous Matsubara 
single-particle Green's functions are obtained at any 



G^{k,uj) 



uj + if] - ^(k) - crii(k. 



-iry) 



(39) 



(77 being a positive infinitesimal) represents the retarded 
{R) single-particle Green's function (for real oj) associ- 
ated with the Matsubara Green's function (|S7|) . since it 
satisfies the requirements of the Baym-Mermin theorem^ 
for the analytic continuation from the Matsubara Green's 
function. 

The first step of the above program is to find the an- 
alytic continuation of S]ii(k, cug) (and S22(k, tOs)) off the 
real axis in the complex z-plane. To this end, it is con- 
venient to express Eii(k, w^) via the spectral form: 



/+00 
-00 



duj' h{k,uj') 



(40) 



With the replacement iuJs — > 2, the spectral representa- 
tion 1)40(1 defines an analytic function Eii(k, z) off the 
real axis. In the case of interest with I]ii(k, los) given by 
Eq. lO, the function hik.uj) of Eq. ||^ reads: 



h{k,oj) = 



dq 
(2^ ^"''1 



{<_klmrfi(q,c. + i?(q-k)) 



X [f{Ei<i-k))+b{u; + E{q-k))] 

+ t.2_kImFfi(q,Lj-S(q-k)) 

X [f{^E{q~k))+b{u;-E{ci-k))]} (41) 
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where f{x) = [exp(j3x) + 1]~^ is the Fermi distribution 
while u^ and w^ are the BCS coherence factors. To obtain 
the expression l|41|l . a spectral representation has been 
also introduced for Fn entering Eq. I|13|) . by writing: 



1 r+°° 

I" J-oo 



,I niFfi(q,c^O 
iQij — uj' 



(42) 



Here, the spectral function F{\(q, a;) is defined by 
Fii(q, iVii, — > w + irf)^ which is obtained from the defini- 
tions 0-101) with the replacement iO,^ -^ uj + irj after the 
sum over the internal frequency a;„ has been performed 
therein. Even in the absence of an explicit Lehmann 
representation for Fn, in fact, it can be shown that the 
spectral representation H42|l holds provided the function 
Fii(q, iili, — > z) of the complex variable z is analytic 
off the real axis. The crucial point is to verify that the 
denominator in Eq. |j2Jl with the replacement ifli, -^ z 
never vanishes off the real axis. This property can be ex- 
plicitly verified in the strong-coupling limit, as discussed 
below. For arbitrary coupling, we have checked it with 
the help of numerical calculations. For the validity of 
the expression l|l^ . it is also required that Fii(q, z) van- 
ishes for \z\ — > cxD. This property can be proved directly 
from Eqs. ©"O, according to which Fii(q, 2;) has the 
asymptotic expression 



rii(q,2:) 



-1 



l3/2 



(43) 



-z + 



2fi 



and thus vanishes for \z\ -^ 00. Once Eii(k, z) has been 
explicitly constructed according to the above prescrip- 
tions, S22(k, z) is obtained as — Eii(k, —z) in accordance 
with Eq. (|T2Jl. 

From the spectral representation H40() for I]ii(k, z), it 
can be further shown that Eii(k, 2) vanishes when l^j — > 
c» along any straight line parallel to the real axis with 
Im z y^ 0. It can also be shown that Iml]ii(k, z) < 
(ImEii(k, z) > 0) when Imz > (Imz < 0). This prop- 
erty follows from the spectral representation of I]ii(k, z), 
provided ft,(k, cj) > in Eq. (|40|l . For arbitrary coupling, 
we have verified that /i(k, w) > with the help of nu- 
merical calculations. In the strong-coupling limit, this 
condition can be explicitly proved, as discussed below. 

From these properties of Sii(k, z) (and S22(k, z)) it 
can then be verified that the function 



crii(k, z) = i;ii(k, z) 



A2 



z + e(k)+ Sn(k,-z) 



(44) 



satisfies the requirements (i)-(iv) stated after Eq. (|38|l . 
With the replacement z — > w -f irj, Eq. H39|) follows even- 
tually on the real frequency axis for the retarded Green's 
function G'^(k, w). 

For later convenience, we introduce the following no- 
tation on the real frequency axis : 



such that Eii(k, — oj — irj) — T,fi(k, —uj)* and 
o-fi(k,Lj) = aii{k,uj + iT]) 
= Efi(k,a;) 



uj + ii] + ^(k)+ S]«(-k,-w) 



(46) 



From Eq. iPHjl it is also clear that ImSf\(k, w) = 
— /i(k, cj), and that ReSf^(k, w) and ImEfj(k, w) are re- 
lated by a Kramers-Kronig transform. 

As anticipated, the properties of the function Eii(k, z), 
required above to obtain the retarded Green's function 
(|39|l on the real axis, can be explicitly verified in the 
strong-coupling limit without recourse to numerical cal- 
culations. In this case, the approximate expression © 
can be used for Fn. This can be cast in the form 142|l . 
with 

ImFfi(q,c^) = -^^[vl{q)S{uj + EB{q)) 

-ul{ci)S{uo-EB{ci))]-m 

Entering the expression (|47|l into Eq. H41(l and the result- 
ing expression into Eq. (|40|) . one obtains for Eii(k, cj^) 
the sum of four terms: 



i;ii(k, Ws) 



dq 



(27r)3 



Stt 
2 . ^ 2 K^i3(q)) + /(i?(q-k)) 



1, and 



Efi(k,t^) = Eii(k,t^-f iT?) 



(45) 



+ "^^'^^'''»-'^SB(q)+i^(q-k)-zc.. 
b{-Es{ci))+f{E{q^\.)) 

+ ''^^'^^"'>-''^MiyT^(^^n^yT7^ 

+ v^(.)v' K-i^B(q)) + /(-g(q-k)) ' 

Since in strong coupling /(_E(k)) -^ 0, u^ 
w^ ^ 0, the second and fourth term within braces on the 
right-hand side of the Matsubara expression H48|) may be 
dropped. The simplified expression (|21() then results from 
Eq. (|48|l . In the strong-coupling limit, one would then be 
tempted to perform the analytic continuation ilUs — + z 
directly from the expression (|21|1 . Care must, however, 
be exerted on this point since the processes of taking the 
strong- coupling limit and performing the analytic contin- 
uation may not commute. By performing the analytic 
continuation iujs -^ z directly in Eq. (|48|l one, in fact, 
obtains two additional terms with respect to the analytic 
continuation of Eq. (|21|l . These two additional terms 
cannot be dropped a priori by the presence of the small 
factor Wp_i(. in the strong-coupling limit, because for real 
z the corresponding energy denominators may vanish. 
Retaining properly these two additional terms indeed af- 
fects in a qualitative way the spectral function A(k, w) 
in the strong-coupling limit, as discussed in Sec. III. 

With the expression obtained by the analytic continu- 
ation iujs — > z of Eq. H48|) , one can prove explicitly that 
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Sii(k, 2:) is analytic off the real axis and vanishes like 
z~^ along any straight line parallel to the real axis with 
Im z ^ 0, and that sgn[IniSii(k, z)] — — sgnpm z]. In 
this way, the properties of the function I]ii(k, z), re- 
quired to obtain the retarded Green's function (|39|l on 
the real axis, are explicitly verified in the strong-coupling 
limit. 

Once the retarded Green's function has been obtained 
in the form (|39|l according to the above prescriptions, its 
imaginary part defines the spectral function 



A(k,Lu) = -(l/7r)ImG'^(k,tj) 



(49) 



which will be calculated numerically in Sec. Ill for a wide 
range of temperatures and couplings. In the Appendix, 
it is shown at a formal level that A(k, oj) satisfies the sum 
rule ljA6p . This sum rule will be considered an important 
test for the numerical calculations of Sec. III. To this 
end, it is necessary to prove that the sum rule ljA6p holds 
even for our approximate theory based on the Dyson's 
equation lfT7|) . 

To prove the sum rule ljA6|l for the approximate theory, 
it is sufficient that the approximate Gn (k, z) (from which 
the retarded Green's function H39|) results when z = uj + 
irf) behaves like z~^ for large \z\. This property is verified 
by our theory, as shown above. As a consequence: 



/-l-OO 1 nOO 

dujA{\<i,u) = Im / <iwG^(k, 

® dajGii(k, z) 



OJ 



-— Im 

TT 



1 (50) 



where the contour G is a half-circle in the upper-half com- 
plex plane with center in the origin, large radius (such 
that the approximation Gii(k, z) ^ z~^ is valid), and 
counterclockwise direction. 

Finally, the analytic continuation of the anomalous 
Matsubara single-particle Green's function H36|l can be 
obtained by following the same procedure adopted for 
the normal Green's function H35|l . One writes for the 
retarded anomalous Green's function 

F«(k,c.)-A[(c. + *r7-e(k)-Sfi(k,c.)) 

X (c. + Z77 + e(k) + Sfi(-k, ^u;)*) - A2]-i . (51) 

in the place of Eqs. (|39|l and H4t)|) . In this case, the an- 
alytic properties of Sii(k, z) {i — 1,2) discussed above 
imply that Gi2(k, z) ~ — A/z^ asymptotically for large 
\z\. As a consequence, the imaginary part of _F^(k, w) 



S(k,w) = -{l/Tr)lmF^ik,Lu) 
satisfies the two following sum rules: 



/ + 00 
dcjB{k,uj) = 
-OO 



and 



+ 00 



div B (k^ cu) Lo = — A 



(52) 



(53) 



(54) 



These sum rules can be verified by introducing the con- 
tour G as in Eq. H5()|l . Note again that these sum 
rules (which are proved on general grounds in the Ap- 
pendix for the exact anomalous retarded single-particle 
Green's function) follow here from our approximate form 
of F^(k,uj) only on the basis of the properties of ana- 
lyticity. Verifying numerically the sum rules H5()|l . H53|l . 
and H54|l at any coupling and temperature will, in prac- 
tice, constitute an important check on the validity of the 
above procedure for the analytic continuation. 

An additional numerical check on the validity of the 
whole procedure at intermediate-to-weak coupling will 
be provided by the merging of the results, obtained by 
calculating the spectral function A(k, cu) when approach- 
ing Tc from below, with the results previously obtained 
in the normal phaseS^ when approaching Tc from above. 



III. NUMERICAL RESULTS AND DISCUSSION 

In this section we present the numerical results based 
on the formal theory developed in Sec. II. Specifically, in 
Sec. IIIA we present the results obtained by solving the 
coupled equations H19|l and H2(J|I for the order parameter 
and chemical potential. Section IIIB deals instead with 
the numerical calculation of the spectral function H49|) 
in the broken-symmetry phase, over the whole coupling 
range from weak to strong. 



A. Order parameter and chemical potential 

Before presenting the numerical results for A and /z, 
it is worth outlining briefly the numerical procedure we 
have adopted. 

At given temperature and coupling, the coupled equa- 
tions H19() and (|20() for the unknowns A and fj, are solved 
via the Newton's method. This requires knowledge of the 
self-energy Sii(k,a;s) of Eq. ((T^ . with rii(g) obtained 
from Eqs. ©-Q. [As anticipated, in the numerical cal- 
culations we neglect T,^2 ^^ comparison to 2^2^^, since 
inclusion of £{"2 is not required to recover the Bogoli- 
ubov results in the strong-coupling limit, as shown in 
Sec. IID.] 

To this end, the frequency sums in Eqs. ^ and Q are 
evaluated analytically, while the remaining wave-vector 
integral is calculated numerically by the Gauss-Legendre 
method. In particular, the radial wave-vector integral 
extending up to infinity is partitioned into an inner and 
an outer region, with the transformation |p| -^ l/|p| 
exploited in the outer region. 

The bosonic frequency sum in Eq. 113|) requires spe- 
cial care, owing to its slow convergence and the lack of 
an intrinsic energy cutoff within our continuum model. 
We have accordingly partitioned this frequency sum into 
three regions, separated by the frequency scales flci and 
ric2 (with < flci < f^ca)- For |ri,y| < ^a, the frequency 
sum is calculated explicitly. For fici < \^u\ < ^C2i the 
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FIG. 3: Order parameter A (in units of e^) vs temperature 
(in units of Tc) for different values of the coupling [kpcLF)'^ ■ 
Results obtained by the inclusion of fluctuations (full lines) 
are compared with mean-fleld results (dashed lines). 



frequency sum is approximated with great accuracy by 
the corresponding numerical integral, owing to the slow 
dependence of Fn on Q,^. Finally, the tail of the fre- 
quency sum for ric2 < l^fl (where the asymptotic ex- 
pression (|^ yields Fn oc (zfi^)"^/^) is evaluated an- 
alytically. Typically, f^ci is taken of the order of the 
largest among the energy scales jw^l, A, |/i|, q^/(2m), and 
k^/(27n); Q,C2 is then taken at least ten times Vlci- It 
turns out that it is most convenient to apply this pro- 
cedure to the frequency sum in Eq. I|13|) after the inte- 
gration over the two angular variables of the wave vector 
q has been performed analytically; the remaining radial 
wave-vector integration is then performed numerically, 
with a cutoff much larger than the wave-vector scales 
|k|, V2mA, and ^/\iJ\. 
Finally, the frequency sum in the particle number 



equation H19() is evaluated by adding and subtracting the 
BCS Green's function Qn on the right-hand side of that 
equation, in order to speed up the numerical convergence. 
Matsubara frequencies are here summed numerically up 
to a cutoff frequency, beyond which the sum is approxi- 
mated by the corresponding numerical integral. The ra- 
dial part of the wave-vector integral in Eq. (|19|l is also 
calculated numerically up to a cutoff scale beyond which 
a power-law decay sets in, so that the contribution from 
the tail can be calculated analytically. 

With the above numerical prescriptions, we have ob- 
tained the behavior of A and /i vs temperature and cou- 
pling reported in Figs. 13161 

Specifically, Fig. O shows the order parameter A 
vs temperature for different couphngs [{kpaF)^^ = 
—0.5,0.5, 1.2, from top to bottom], in the window —1 < 
{kpaF)^^ ^ +1 where the crossover from weak to strong 
coupling is exhausted. Comparison is made with the cor- 
responding curves obtained within mean field (dashed 
lines), when the BCS Green's function Qn enters Eq. ifliHl 
in the place of the dressed Gn. In these plots, the 
temperature is normalized with respect to the critical 
temperature T^ for the given coupling. This compari- 
son shows that fluctuation corrections on top of mean 
field get progressively important at given coupling as the 
temperature is raised toward T,.. Close to Tc, fluctuation 
corrections become even more important upon approach- 
ing the strong-coupling limit. Near zero temperature, on 
the other hand, fluctuation corrections become negligi- 
ble when approaching strong coupling. This confirms the 
expectation that, near zero temperature, the BCS mean 
field should be rather accurate both in the weak- and 
strong-coupling limits^ . 

Note from Fig.Othat A jumps discontinuously close to 
the critical temperature when fluctuations are included 
on top of mean field. This jump becomes more evident 
as the coupling is increased. It reflects an analogous be- 
havior of the condensate density near the critical tem- 
perature as obtained by the Bogoliubov theory for point- 
like bosons^ . In the present theory this jump is carried 
over to the composite bosons, even at fermionic couplings 
[as in the middle panel of Fig. O ] when the composite 
bosons are not yet fully developed. When the fermionic 
coupling increases beyond the values reported in Fig. O 
however, the residual interaction between the composite 
bosons decreases further and the jump becomes progres- 
sively smaller. More reflned theories for point-like bosons 
(see, e.g., Ref. WS) remove the jump of the bosonic con- 
densate density, which thus should be considered as an 
artifact of the Bogoliubov approximation. Apart from 
this jump, note that when the temperature is decreased 
below Tc the order parameter A grows more rapidly with 
the inclusion of fluctuations than within mean field. 

Figure ^ shows the chemical potential // vs tempera- 
ture for the same coupling values of Fig. |3| Note that in 
weak coupling the chemical potential decreases slightly 
upon moving deep in the superconducting phase from 
Tc to r = 0, in agreement with the BCS behavior. In 
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FIG. 4: Chemical potential /i (in units of ep) vs temperature 
(in units of Tc), for the same values of the coupling {kpciF)'^ 
as in Fig. 13 
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FIG. 5: Order parameter A at T = (in units of ep) vs 
the coupling {kpap)' . Results obtained by the inclusion of 
fluctuations (full line) are compared with mean-field results 
(dashed line). 



FIG. 6: Chemical potential /i at T = (in units of ep for 
/i > and of eo/2 for /i < 0) vs the coupling {kpap)~^ . 
Results obtained by the inclusion of fluctuations (full line) 
are compared with mean-field results (dashed line). 



at zero temperature. 

Note, finally, that the values for A and /i obtained 
from our theory at T = with the coupling value 
(kpCLp)^^ — are in remarkable agreement with a recent 
Quantum Monte Carlo calculation'*'^ performed for the 
same coupling. Our calculation yields, in fact, A/ep — 
0.53 and fi/ep = 0.445, to be compared with the val- 
ues A/ep = 0.54 and n/ep = 0.44 ± 0.01 of Ref. IH 
[In contrast, BCS mean field yields A/ep = 0.69 and 
fi/ep = 0.59.] 

In summary, the above results have shown that, for 
thermodynamic quantities like A and /z, fluctuation cor- 
rections to mean-field values in the broken-symmetry 
phase are important only as far as the temperature de- 
pendence is concerned, while at zero temperature the 
mean-field results are reliable. 



strong coupling the opposite occurs, reflecting the behav- 
ior of the bosonic chemical potential /is — 2/i -I- eg within 
the Bogoliubov theory. It should be, however, men- 
tioned that with improved bosonic approximations^^, the 
bosonic chemical potential would rather decrease upon 
entering the condensed phase. 

Figure shows the order parameter A at zero temper- 
ature (full line) and the corresponding mean-field value 
(dashed line) vs the coupling {kFO,F)~^- While A in- 
creases monotonically in absolute value from weak to 
strong coupling (as expected on physical grounds), the 
relative importance of the fiuctuation corrections to the 
order parameter at zero temperature (over and above 
mean field) reaches a maximum in the intermediate- 
coupling region, never exceeding about 30%. This results 
confirms again that the BCS mean field is a reasonable 
approximation to the ground state for all couplings. 

Figure 1^ shows the chemical potential fi at zero tem- 
perature vs the coupling parameter (kpaF)"^- The re- 
sults obtained by the inclusion of fluctuations (full lines) 
are compared with mean field (dashed lines). Even for 
this thermodynamic quantity the fluctuation corrections 
to the mean-field results appear to be not too important 



B. Spectral function 

For a generic value of the coupling, calculation of the 
imaginary part of the retarded self-energy ImS{\ (k, oj) = 
— /i(k, w) (with /i(k,w) given by Eq. H41() ) requires us to 
obtain the imaginary part of the particle-particle lad- 
der r{\ on the real- frequency axis, as determined by the 
formal replacement i^i, —^uj + irjin the Matsubara ex- 
pressions (01-0). After performing the frequency sum 
therein, the wave-vector integrals of Eqs. Q and l@J for 
the functions Xy (Qj ^^i' ~^ i^+ir]) iihj) = Ij 2) are evalu- 
ated numerically, by exploiting the properties of the delta 
function for the imaginary part and keeping a finite albeit 
small value of 77 {— lO^^y/ji'^ + A^) for the real part. 

Direct numerical calculation of the imaginary part of 
the particle-particle ladder fails, however, when this part 
has the structure of a delta function for real u at given 
q. This occurs when the determinant in the denomi- 
nator of Eq. ^ vanishes for real cu. To deal with this 
delta function, let's first consider the case T = for 
which three cases can be distinguished, according to: (i) 
ap < and fi > (weak-to- intermediate coupling); (ii) 
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FIG. 7: Dispersion c^(q) of the pole of rfi(q,(.j) at T = 
(full lines) and boundary of the particle-particle continuum 
(dashed lines) for three characteristic couplings. 



FIG. 8; Imaginary part of the self-energy a^^ for |kj — k^i 
vs frequency (in units of ep) at different temperatures for the 
coupling values (kpaF)"^ = -0.5 (a), 0.1 (b), and 0.5 (c). 



ap > and /i > (intermediate coupling); (iii) ap > 
and n < (intermediate-to-strong coupling). The curves 
ix'(q) where the (analytic continuation of the) determi- 
nant in the denominator of Eq. ^ vanishes are shown 
(full hues) for these three cases in Figs. [7|(a), (b), and 
(c) , respectively. In these figures we also show the bound- 
aries (dashed lines) delimiting the particle-particle con- 
tinuum, where the imaginary part of the particle-particle 
ladder is nonvanishing and regular (in the sense that it 
does not have the structure of a delta function) . At finite 
temperature, the sharp boundary of the particle-particle 
continuum smears out, owing to the presence of Fermi 
functions after performing the sum over the Matsubara 
frequencies in Eqs. Q and @. The Fermi functions pro- 
duce, in fact, a finite (albeit exponentially small with 
temperature) imaginary part of the particle-particle lad- 
der also below the (dashed) boundaries of Fig.[71 resulting 



in a Landau-type damping of the Bogoliubov- Anderson 
mode (^(q)- In addition, the finite imaginary part broad- 
ens the delta-function structure centered at the curves 
Ci;(q) of Fig. {7\ turning it into a Lorentzian function. 
In practice, our numerical calculation takes advantage 
of this broadening occurring at finite temperature, and 
deals with smooth Lorentzian functions instead of the 
delta-function peaks.— 

As a further consistency check on our numerical cal- 
culations, we have sistematically verified that the three 
sum rules (I5UII . H53|l . (|54|) are satisfied within numeri- 
cal accuracy, for all temperatures and couplings we have 
considered. 

The imaginary and real parts of the retarded self- 
energy a^i (k, Lu) obtained from Eq. (|45|l are shown, re- 
spectively, in Figs. IHlandOas functions of frequency at 
different temperatures and for different couplings (about 
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FIG. 9: Real part of the self-energy (7j\ for |k| — fc^/ vs 
frequency (in units of ep) at different temperatures for the 
coupling values (kpaF)^^ ~ -0.5 (a), 0.1 (b), and 0.5 (c). 



FIG. 10; Spectral function for |k| — k^^i vs frequency (in 
units of ef) at different temperatures for the coupling values 
[kFap)-^ = -0.5 (a), 0.1 (b), and 0.5 (c). 



the crossover region of interest). The magnitude of the 
wave vector k is taken in Figs.|HlandOat a special value 
(denoted by fc^'), which is identified from the behav- 
ior of the ensuing spectral function A{]s.,ui) when per- 
forming a scanning over the wave vector (see Fig. 12 
below). Accordingly, fc^/ is chosen to minimize the 
gap in the spectral function, in agreement with a stan- 
dard procedure in the ARPES literature. On the weak- 
coupling side (when the the self-energy shift Sq discussed 
in Sec. IIB is included in our calculation), k^i coincides 
with ^J2m{ii — Eg). On the strong-coupling side (when 
H becomes negative) one takes instead fc^' =0. 

For all couplings here considered, the progressive evo- 
lution found in j4(k, uj) (from the presence of a pseudogap 
about a; = at Tc to the occurrence of a superconducting 
gap near zero temperature) stems from the interplay of 



the two contributions in Eq. (|46|l to the imaginary part 
of crfi (k,a;) about lo — 0. Specifically, for intermediate- 
to-weak coupling (with /i > 0) the first term on the right- 
hand side of Eq. 1)46(1 (which is responsible for the pseudo- 
gap suppression in j4(k, uj) at Tc) would produce a narrow 
peak structure in A{k^i , w) about w = upon lowering T, 
since ReE^]^(fc^', w) — Sq vanishes while |ImE(\(fc^',Ci;)| 
becomes progressively smaller. The presence of the sec- 
ond term on the right-hand-side of Eq. (|46|l . however, 
gives rise to a narrow peak in Im cr{\ (fc^/ , lu) about tj = 0, 
as seen from Fig. |S1 (a), resulting in a depression of 
A(k^i,ijj) about lo — 0. [This occurs barring a small 
temperature range close to Tc^ where the second term on 
the right-hand-side of Eq. (|^ is not yet well developed.] 
At larger couplings (when ^ < 0), the first term on the 
right-hand-side of Eq. (|46|l would not produce a peak in 
yl(k = 0, w) about w = upon lowering the temperature. 
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FIG. 11: Temperature dependence of the total weight of the 
spectral function at negative frequencies for different coupling 
values about the crossover region. 



because |/i| +ReSf;^(fc = O^u) does not correspondingly 
vanish in this case even though Im S{\ (k = 0, oj) does. In 
addition, in this case the second term on the right-hand- 
side of Eq. (|46() does not produce a peak in yl(k ~ 0, uj) 
about a; = 0. 

FigurelTUIshows the resuhing spectral function A(k, oj) 
vs Lu for |k| = k^i at different temperatures and cou- 
plings. In all cases, at Tc there occurs only a broad 
pseudogap feature both for uj > and to < 0. [For 
photoemission experiments only the case cj < is rel- 
evant, so that we shall mostly comment on this case in 
the following.] A coherent peak is seen to grow on top 
of this broad pseudogap feature as the temperature is 
lowered below Tc. When zero temperature is eventually 
reached, the pseudogap feature is partially suppressed in 
favor of the coherent peak, which thus absorbs a substan- 
tial portion of the spectral intensity. This interplay be- 
tween the broad pseudogap feature and the sharp coher- 
ent peak results in a characteristic peak- dip-hump struc- 
ture, which is best recognized from the features for weak- 
to-intermediate coupling. Generally speaking, this coher- 
ent peak (and its corresponding counterpart at positive 
frequencies) for intermediate-to-weak coupling is associ- 
ated with the two dips in Im crfj^ symmetrically located 
about zero frequency [cf. Figs.|Hl(a) and (b)]. In strong 
coupling, instead, the coherent peak results from a deli- 
cate balance between the real and imaginary parts of cr(]^ 
near the boundary of the region where ImE{\ — 0. 

An interesting fact is that the weights of the nega- 
tive and positive frequency parts of the spectrum turn 
out to be separately (albeit approximatively) constant 
as functions of temperature for given coupling, as shown 
in Fig. E] for three characteristic couplings. This implies 
that, for a given coupling, the coherent peak for tj < 
grows at the expenses of the accompanying broad pseu- 
dogap feature upon decreasing the temperature. 

The result that the total area for negative uj should 
be (approximately) constant as a function of tempera- 
ture can be realized also from the analytic results in the 
extreme strong-coupling limit discussed in Sec. IID. Tak- 
ing the analytic continuation of the Matsubara Green's 



function H29() (which is appropriate in the strong-coupling 
limit as far as this total area is concerned, as it will be 
shown below) results, in fact, in the total weight f;^(k) 
of the UJ < Q region being independent of temperature, 
since the combination A^ -I- Aq entering the expression 
of i^(k) is proportional to the total density in this limit 
[cf. Eq. 132)] ■ 

Returning to Fig. 1101 it is also interesting to comment 
on the positions of the pseudogap feature and the coher- 
ent peak as functions of temperature for given coupling. 
The position of the coherent peak depends markedly on 
temperature, shifting progressively toward more negative 
frequencies as the temperature is lowered. In particu- 
lar, for weak-to-intermediate coupling the position of the 
coherent peak about coincides with (minus) the value 
of the order parameter A. In the strong-coupling re- 
gion (where /i < 0), on the other hand, its position is 
about at — yA^ + zi^- This remark entails the possi- 
bility of extracting two important quantities from the 
temperature evolution of the coherent peak in the spec- 
tral function: (i) The frequency position of this peak 
when approaching T^ determines whether /i is positive 
(when the peak position approaches uj — Q) oy negative 
(when the peak position approaches — ImI), corresponding 
to weak-to-intermediate coupling and strong coupling, 
respectively; (ii) In both cases, the temperature depen- 
dence of the order parameter can be extracted from the 
frequency position of the coherent peak. 

The above results for the coherent peak contrast some- 
what with the position of the pseudogap feature by de- 
creasing temperature below Tc, also determined from 
Fig. ^1 The broad pseudogap feature does not depend 
sensitively on temperature for all couplings shown in this 
figure. This indicates that the broad pseudogap feature 
does not relate to the order parameter below Tc- 

As far as the spectral function is concerned, one of 
the key results of our theory is thus the presence of two 
structures (coherent peak and pseudogap) , which behave 
rather independently from each other as functions of tem- 
perature and coupling. This result, which is also evi- 
denced by the behavior of the experimental spectra in 
tunneling experiments on cuprates^, originates in our 
theory from the presence of two distinct contributions to 
the self-energy, namely, the BCS and fluctuation contri- 
butions of Eq. (fT7|l . While the broad pseudogap feature 
at T < Tc develops with continuity from the only feature 
present at T > Tc, the coherent peak per se would be 
present in a BCS approach even in the absence of the 
fluctuation contribution. This remark, of course, does 
not imply that the two contributions to the self-energy 
of Eq. (|17|l are totally independent from each other. They 
both depend, in fact, on the value of the order parame- 
ter A which is, in turn, determined by both self-energy 
contributions via the chemical potential. 

A further important feature that can be extracted from 
our calculation of the spectral function is the evolution of 
the coherent peak for varying wave vector at fixed tem- 
perature and coupling. Figure [T^ reports A(k,uj) vs uj 
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FIG. 12: Spectral function at different wave vectors about fc^j/ 
for T = O.GTc vs frequency (in units of ep) for the coupling 
value {kFaF)~ = —0.5. 



for different values of the ratio k/kf^i about unity when 
{kpaFy^ = -0.5 and T/Tc = 0.6. Here, fc/fc^, = 1 
identifies the underlying Fermi surface that represents 
the "locus of minimum gap". When fc/fc^/ < 1, there is 
a strong asymmetry between the two coherent peaks at 
Lu < and w > 0, with the peak at cj < absorbing 
most of the total weight. The situation is reversed when 
k/kfj,' > 1. When fc/fc^' — 1 the spectrum is (about) 
symmetric between lu and —lu. In addition, when follow- 
ing the position of the coherent peak at w < starting 
from k/k^i < 1, one sees that this position moves toward 
increasing lu, reaches a minimum distance from lu — 0, 
and bounces eventually back to more negative values of 
LU. The value of the minimum distance from lu — iden- 
tifies an energy scale A^. At the same time, the weight 
of the coherent peak at w < progressively decreases for 
increasing fc/fc^/ starting from k/k^i < 1. When k/k^i 
becomes larger than unity, the weight of the coherent 
peak is transferred from negative to positive frequencies. 
This situation is characteristic of the BCS theory, where 
only the coherent peaks are present without the accompa- 
nying broad pseudogap features. Our calculation shows 
that this situation persists also for couplings values in- 
side the crossover region, where the presence of the pseu- 
dogap feature is well manifest due to strong supercon- 
ducting fluctuations. [Sufficiently far from the underly- 
ing Fermi surface, the coherent peak and the pseudogap 
feature merge into a single structure, as it is evident from 
Fig. El In this case, the above as well as the following 
considerations apply to the structure as a whole and not 
to its individual components.] 

Figure lTST a) summarizes this finding for the dispersion 
of the coherent peaks, by showing the positions of the 
two coherent peaks as extracted from Fig. E]vs k/k^'. 
These positions are compared with the two branches 
±A/^(k)^ -t- A^ of a BCS-like dispersion, where A.^ is 
also identified from Fig. ^1 [The value of A^ turns out 
to about coincide with the value of the order parameter 
A at the same temperature, see below.] The correspond- 
ing evolution of the weights of these peaks is shown in 
Fig. ll^r b). where the characteristic feature of an avoided 
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FIG. 13: (a) Positions of the coherent peaks (in units of e^) vs 
the wave vector as extracted from Fig. 1121 . Positive (squares) 
and negative (circles) branches are compared with BCS-like 
dispersions (full lines), as explained in the text, (b) Corre- 
sponding weights vs the wave vector, with particle-like (full 
line) and hole-like (dashed line) contributions. 



crossing is evidenced. The dispersion of the positions and 
weights of the coherent peaks shown in Fig. 1 131 compare 
favorably with those recently obtained experimentallyi& 
for slightly overdoped Bi2223 samples below the critical 
temperature (for T/Tc — 0.6). 

An additional outcome of our calculation is reported 
in Fig. El where the distance A^ of the coherent peak 
in A{\<i,uj) from a; = at |k| = fc^' is compared at low 
temperature with the order parameter A when /z > and 
with y^i^ + A^ when /i < 0. This plot thus compares dy- 
namical and thermodynamic quantities. The good agree- 
ment between the two curves confirms our identification 
of the coherent-peak position in ^(k, w) with the mini- 
mum value of the excitations in the single-particle spec- 
tra according to a BCS-like expression (where the value 
of the order parameter A is, however, obtained by in- 
cluding also fiuctuation contributions). 

Finally, it is interesting to comment on the strong- 
coupling result (|28|l for the diagonal Green's function, 
with a characteristic double-fraction structure. The cor- 
responding spectral function A(k, uj), obtained from that 
expression after performing the analytic continuation 
iojn — > UJ -\- irj, shows only a single feature for w < 0, 
with a temperature-independent position. This contrasts 
the numerical results we have presented [cf. in partic- 
ular Fig. El- This difference is due to the fact that. 
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FIG. 14: Position Am (in units of tp) of the quasi-particle 
peak at T = O.lTc vs the coupling (/cFflF)"^ (full line). The 
dashed line corresponds to the value of the order parameter 
A when /i > and of ^A^ + p? when /i < 0. 



FIG. 15: Spectral function vs frequency for (fcFflF)"^ ~ 2.0 
and T jTc = 0.1, obtained by taking alternatively the analytic 
continuation of En from the expression l|48|l (full line) or from 
the expression 1231 (broken line). 



in our numerical calculation, the analytic continuation 
has been properly performed before taking the strong- 
coupling limit, as emphasized in Sec. HE. With this pro- 
cedure, in fact, the pseudogap structure and the coherent 
peak remain distinct from each other even in the strong- 
coupling limit, without getting lumped into a single fea- 
ture. Such a noncommutativity of the processes of taking 
the analytic continuation and the strong-coupling limit 
was noted already in a previous paper^^ when studying 
the spectral function above Tc- More generally, the oc- 
currence of this noncommutativity is expected whenever 
one considers approximate expressions in the Matsub- 
ara representation and takes the analytic continuation of 
these expressions to real frequency. 

To make evident the noncommutativity of the two pro- 
cesses, we show in Fig. ^1 the spectral function yl(k = 
0,w) for (kFap)'^ = 2.0 and T/Tc = 0.1, obtained 
by two alternative methods, namely: (i) Using the an- 
alytic continuation of the expression H48|) for Sn where 
iu>s —> Lo + irj (full line); (ii) Taking the strong-coupling 
expression (|23|l for Sn, in which the analytic continua- 
tion iLUg —I- Lu + irj is performed (broken line). Method 
(i) results in the presence of two distinct structures in 
yl(k, w) for cu < 0, corresponding to the coherent (delta- 
like) peak and the broad pseudogap feature. Method (ii) 
gives instead a single delta-like peak. It is interesting 
to note that the total spectral weight of the two peaks 
for w < obtained by method (i) (=0.049 for the cou- 
pling of Fig. I15|l about coincides with the weight of the 
delta-like peak (=0.044) obtained by method (ii). [We 
have verified that this correspondence between the spec- 
tral weights persists also at stronger couplings.] 

These remarks explain the occurrence of a single fea- 
ture in the spectral function as obtained by a different 
theory based on a preformed-pair scenario^. In that 
theory, a single-particle Green's function with a double- 
fraction structure is considered in the Matsubara repre- 
sentation for any coupling, and correspondingly a single 
feature in the spectral function is obtained for real fre- 
quencies^. Our theory shows instead the appearance of 



two distinct energy scales (pseudogap and order param- 
eter) in the spectral function below Tc- 

We are thus led to conclude that the occurrence of 
two distinct energy scales below Tc in photoemission and 
tunneling spectra should not be necessarily associated 
with the presence of an "extrinsic" pseudogap due to ad- 
ditional non-pairing mechanisms, as sometimes reported 
in the literature^. 



IV. CONCLUDING REMARKS 

In this paper, we have extended the study of the BCS- 
BEC crossover to finite temperatures below Tc. This has 
required us to include (pairing) fluctuation effects in the 
broken-symmetry phase on top of mean field. Our ap- 
proximations have been conceived to describe both a sys- 
tem of superconducting fermions in weak coupling and a 
system of condensed composite bosons in strong coupling, 
via the simplest theoretical approaches valid in the two 
limits. These are the BCS mean field (plus superconduct- 
ing fluctuations) in weak coupling and the Bogoliubov 
approximation in strong coupling. To this end, analytic 
results have been specifically obtained in strong coupling 
from our general expression of the fermionic self-energy. 

Results of numerical calculations have been presented 
both for thermodynamic and dynamical quantities. The 
latter have been defined by a careful analytic continua- 
tion in the frequency domain. In this context, a noncom- 
mutativity of the analytic continuation and the strong- 
coupling limit has been pointed out. 

Results for thermodynamic quantities (like the order 
parameter and chemical potential) have shown that the 
effects of pairing fluctuations over and above the BCS 
mean field become essentially irrelevant in the zero- 
temperature limit, even in strong coupling. Results for 
a dynamical quantity like A(k, w) have shown, in ad- 
dition, that two structures (a broad pseudogap feature 
that survives above Tc and a strong coherent peak which 
emerges only below Tc) are present simultaneously, and 
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that their temperature and couphng behaviors are rather 
(even though not completely) independent from each 
other. 

These features produced in the spectral function by our 
theory originate from a totally intrinsic effect, namely, 
the occurrence of a strong attractive interaction (irre- 
spective of its origin). Additional features produced by 
other extrinsic effects could obviously be added on top 
of the intrinsic effects here considered. 

Similar results have recently been obtained in Ref. l49l 
using a boson-fcrmion model for precursor pairing below 
Tc- In that reference, a two-peak structure for A(k, a;) 
has also been obtained, although with a self-energy cor- 
rection introduced by a totally different method. 

The attractive interaction adopted in this paper is the 
simplest one that can be considered, depending on a sin- 
gle parameter only. Detailed comparison of the results 
of this theory with experiments on cuprates would then 
require one to specify the dependence of this effective 
parameter on temperature and doping. 

The simplified model that we have adopted in this 
paper should instead be considered realistic enough for 
studying theoretically the BCS-BEC crossover for Fermi 
atoms in a trap. The occurrence of this crossover in these 
systems is being rather actively studied experimentally at 
presentii^ In this case, the calculation should also take 
into account the external trapping potential by consider- 
ing, e.g., a local version of our theory with local values 
of the density and chemical potential in the trap:^. 
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APPENDIX A: ANALYTIC CONTINUATION 

FOR THE FERMIONIC RETARDED 

SINGLE-PARTICLE GREEN'S FUNCTIONS AND 

SUM RULES BELOW THE CRITICAL 

TEMPERATURE 

In this appendix, we extend below the critical tem- 
perature a standard procedure for obtaining at a for- 
mal level the fermionic retarded single-particle Green's 
functions via analytic continuation from their Matsub- 
ara counterparts. This is done in terms of the Lehmann 
representation^^ and of the Baym-Mermin theoremiS. In 
this context, besides the usual sum rule that holds also 
above the critical temperature^*^, we will obtain two ad- 
ditional sum rules that hold specifically below the critical 
temperature. 

The results proved in this appendix hold exactly, ir- 
respective of the approximations adopted for the Mat- 
subara self-energy. To satisfy the above three sum rules 
with an approximate choice of the self-energy, however. 



it is not required for the ensuing approximation to the 
fermionic single-particle Green's functions to be "con- 
serving" in the Baym's sensed. Rather, it is sufficient 
that the analytic continuation from the Matsubara fre- 
quencies to the real frequency axis is taken properly, as 
demonstrated in Sec. HE with the specific choice Hl(i(l of 
the self-energy. 

We begin by considering the fermionic "normal" and 
"anomalous" retarded single-particle Green's functions in 
the broken-symmetry phase, defined respectively by 

G«(r,t;r',0 = -*^(t-i')({V'T(r,t),^|(r',t')}>(Al) 
F«(r, i; r', t') = ~t0{t - t'){{ij^{r, t),^^', t')}) . (A2) 

Here, d{t) is the unit step function, ijjcr{r,t) is the 
fermionic field operator with spin a ~ (Tii) ^-t po- 
sition r and (real) time t (such that tpa-i^jt) = 
exp{iKt)tl;a{j^) exp{—iKt) with K = H — fj.N in terms 
of the system Hamiltonian H and the particle number 
iV), the braces represent an anticommutator, and (• • •) 
stands for the grand-canonical thermal average. 

The Matsubara counterparts of ljAl|l and (|A2|) are sim- 
ilarly defined by 

G(r,r;r',T') = -{T,[^/j^{r,T)4^\{r' ,r')\) (A3) 
F(r,T;r',r') = -{Tr[^^{r,r)^^{r\r')]), (A4) 

where now ipa-(j,T) = exp(_?CT)'0(j(r) exp(— iiTr), 
'(/'].(r, r) = exp(i4rT)'0j.(r) exp(— iiTr), and Tt is the time- 
ordering operator for imaginary time r. 

The Lehmann analysis for the normal function G^ in 
the broken-symmetry phase proceeds along similar lines 
as for the normal phased. The result is that (for a ho- 
mogeneous system) the wave-vector and (real) frequency 
Fourier transform can be obtained by the spectral repre- 
sentation 



G^(k,uj) 



"°°..'. ^("^'-'^ 



LU — Lu + ir] 



(A5) 



r] being a positive infinitesimal. Here, the real 
and positive definite spectral function A(k, uj) = 
— (l/7r)ImG-^(k, w) satisfies the sum rule 



-j-oc 



dujA(k,uj) = 1 



(A6) 



for any given k, as a consequence of the canonical anti- 
commutation relation of the field operators. 

A similar analysis for the Matsubara normal Green's 
function leads to the spectral representation 



G(k,Ws) = Gii(k,ujs 



^-^'Ai^, (A7) 

— oo 



lOJs 



to 



in terms of the same spectral function A(k, uj) of 
Eq. IjASp . where Us = {2s + l)n/(3 (s integer) is a 
fermionic Matsubara frequency and the diagonal Nambu 
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Green's function has been introduced. The spectral rep- 
resentations IjASfl and (|A7() . together with knowledge of 
the asymptotic behavior G^(k,u!) ~ uj^^ for large |w|, 
are sufficient to guarantee that the retarded normal func- 
tion is the correct analytic continuation of its Matsub- 
ara counterpart in the upper-half of the complex fre- 
quency planed, in accordance with the Baym-Mermin 
theorem^. 

The above Lehmann analysis can be extended to the 
anomalous function (|A2p as well. One obtains 



F^(k,uj) 



duj' ■ 



S(k,c.') 



I T] 



(A8) 



in the place of Eq. IjASp . The new spectral function 
-B(k, Lu) vanishes for large |w| but, in general, is no longer 
real and positive definite. [One obtains for i3(k, w) the 
same formal expression^?, for A(k,uj) in terms of the 
eigenstates \n > of the operators H and N, apart from 
the replacement of |(n'|'(/'|(r = 0)|n)p with {n\tpi{r = 
0)\n') (n'|'(/'|(r = 0)\n).]^ It can then be readily verified 
that B{'k,uj) satisfies the sum rule 



dLuB(k,uj) = 



(A9) 



which is again a consequence of the canonical anticom- 
mutation relation of the field operators. The above prop- 
erties guarantee that F^CkjUj) vanishes faster than lu~^ 
for large \uj\. 

By a similar token, considering the Matsubara anoma- 
lous Green's function leads to the spectral representation 



/ + 00 
duj' 
-oo 



, B{k,u'] 



lUJs 



7, (AlO) 



where the off-diagonal Nambu Green's function has been 
introduced. These considerations suffice again to guar- 
antee that the retarded anomalous function is the correct 
analytic continuation of its Matsubara counterpart in the 
upper-half complex frequency plane, in accordance with 
the Baym-Mermin theorem^. 

Finally, an additional sum rule for _B(p,a;) can be ob- 
tained by using the relation 

dujB{k,uj)uj^i / — F^(k,cj)we-*'^'' 

-oo J-oo 27r 



i dre 



dt 



MO) ) (All) 



and exploiting the equation of motion for the field op- 
erator. For the contact potential we are considering 
throughout this paper, we write 



il '^^'Hr^ ' :Mo) 



= 0^ 



dt 



= -'(;o'5(r)(V'T(r)V'i(r)) 
= -S{r)A (A12) 



in terms of the order parameter A. The expression (|A11|) 
thus becomes: 



/ + 00 
du! B(k, Lj) CO 
-oo 



= -A . 



(A13) 



This constitutes a third sum rule for the spectral func- 
tions in the broken-symmetry phase. 
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